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I consider the problem of defining canonical coordinates and momenta in pure Yang-Mills theory, 
under the condition that Gauss' law is identically satisifed. This involves among other things 
particular boundary conditions for certain dependent variables. These boundary conditions 
are not postulated a priori, but arise as consistency conditions related to the equations of 
motion. It is shown that the theory indeed has a canonical structure, provided one uses a 
special gauge condition, which is a natural generalisation to Yang-Mills theory of the Coulomb 
gauge condition in electrodynamics. The canonical variables and Hamiltonian are explicitly 
constructed. Quantisation of the theory is briefly discussed. 
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1 Introduction 



In this paper I present a novel solution to an old problem in Yang-Mills theory [[j] , namely the 
problem of formulating the Yang-Mills equations of motion as a set of canonical Hamiltonian 
equations. The canonical formulation of Yang-Mills theory has been considered as a necces- 
sary prerequisite for the development of a consistent quantisation procedure of this theory, as 
emphasised e.g. in the textbook on gauge theory by Faddeev and Slavnov 0. I refer to this 
textbook for an account of the state of the art in Yang-Mills theory around 1980. An account 
of more recent developments can be found e.g. in the textbook by Weinberg ||. 

The early attempts to quantise Yang-Mills theory ||, 0] were based on the use of the 
Coulomb gauge condition, familiar from electrodynamics, in Yang-Mills theory. However, as 
shown by Gribov || in 1977, the Coulomb gauge is not a proper gauge condition in Yang-Mills 
theory. It should be noted here that the problems pointed out by Gribov concerning the use of 
the Coulomb condition, the so-called Gribov ambiguities, do not really matter in perturbative 
calculations. 

Another possible approach to a canonical formalism in Yang-Mills theory J/J, is based on 
the so-called Weyl gauge ||, A = 0. However, in this formulation one simply disregards the 
Gauss law, and considers only the canonical Hamiltonian formulation of the remaining Yang- 
Mills equations, which is a rather simple problem. Gauss' law is then in the quantised version 
of this theory introduced by hand, as a condition on the states in the theory. While this may 
be an acceptable approach to quantum Yang-Mills theory, it has, to the best of my knowledge, 
yielded only partial insight into the structure of pure Yang-Mills theory. 

There has recently been an upsurge in interest in non-perturbative aspects of Yang-Mills 
theory, as witnessed e.g. by a paper by Faddeev and Niemi ||. It is emphasised by Faddeev and 
Niemi in their paper, that we are still lacking theoretical understanding of low-energy Yang- 
Mills theory, i.e. the non-perturbative regime, which is relevant for studying questions such as 
colour confinement. The problem of confinement is definitely a question which calls for a proper 
quantum Yang-Mills theory. A proper canonical formulation of the (semiclassical) Yang-Mills 
theory may well be a necessary prerequisite for this, or at least an appropriate starting point. 

The real confinement problem concerns the situation when quark fields are included in the 
formalism, but this extension is not an essential difficulty in the formulation presented here. 

The canonical formulation developed in this paper differs from the Weyl gauge formulation 
in that Gauss' law is implemented identically (in principle). This has necessitated the intro- 
duction of a new gauge condition ||10|| , which is a straightforward and natural generalisation to 
Yang-Mills theory of the ordinary Coulomb gauge condition in electrodynamics. The result- 
ing canonical Yang- Mills theory, in its Schrodinger quantised form, resembles the corresponding 
theory based on the Weyl gauge in many ways, but is apparently not isomorphic to that theory. 
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2 Notation and Conventions 



The basic variables in Yang-Mills theory JTJ are the gauge field G^ u and gauge potential A^, 
respectively. These quantities take values in a convenient (Hermitian) matrix representation 
of the Lie algebra of the gauge group G. The Lie algebra is defined by the structure constants 
f ab ° in the commutator algebra of the Hermitian matrix representatives T a of the Lie algebra 
generators, 

[T a ,T b ] = if ab c T c , (1) 

where appropriate summation over repeated Lie algebra indices (a, b, c, ...,/) is understood. We 
assume G to be semisimple and compact. Then a positive definite Lie algebra metric can be 
given in terms of the following Killing form (h ab ), 

hab — —faV fbd (2) 

The inner product (A, B) of any two Lie algebra valued quantites A = A a T a and B = B a T a is 
defined as follows, 

(A,B) = h ab A a B b (3) 
The form h a b and its inverse h ab are used to lower and raise Lie algebra indices, respectively. 
In the notation introduced so far, we write the gauge potential A^ as follows, 

A,{x) = A;(x)T a (4) 

where the argument x stands for a space-time point in Minkowski space. The gauge field G^ U (A) 
is then given as follows, 

G^(A) = dvA^x) - d„A v {x) - ig[A^x) , A v {x)\, (5) 

where g is an arbitrary nonvanishing real parameter, which is introduced for convenience. An 
alternative to the notation (^) is the following, 

G%(A) = d v Al{x) - d,A a u (x) + gf bc a A%x), A c v {x) (6) 



The remaining notation is fairly conventional. Greek letters fi, ... are spacetime indices 
which take values in the range (0, 1, 2, 3). These indices are lowered (raised) with the standard 
diagonal Minkowski metric g^ u (g^) with signature (+, — , — , — ). Latin indices from the middle 
of the alphabet (A;, £, ...) are used as space indices in the range (1, 2, 3). Unless otherwise stated, 
repeated indices are always summed over, be they Lie algebra-, spacetime- or space indices. 

I finally state the gauge transformation formulae for the gauge potential and field, respec- 
tively, in my notation. Let Q(x) denote a general gauge transformation, the generic form of 
which is a unitary matrix of the following form, 

Q(x) =exp(iga a {x)T a ), (7) 
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where the functions a a (x) are any sufficiently smooth real- valued functions on Minkowski space. 
The relevant gauge transformation formulae are as follows, 

Ap{x) A'^x) = n-\x)A tl {x)Q(x) + l -(d,n~\x)Mx), (8) 

and 

Gfj,v{A) G^{A') = n~\x)G lxv {A)n{x) (9) 



3 The equations of motion and Gauss' law 

The Yang-Mills action S, with a general potential A, which at this stage is not assumed to 
satisfy any particular gauge condition, is as follows, 

S = ~\ J d A x{G, v {A),G^{A)) (10) 



It is convenient to write the action (|T0|) in terms of a Lagrangian L, although this treats time- 
and space indices on an unequal footing. The action S is the integral of a Lagrangian L in an 
appropriate time interval Thus, 



/ dx°L, 

Jx° 



S I <I., V 'L. (II) 
where 

L = -\j v d 3 x (G ok (A), G ok (A)) - - A J v d 3 x (G M (A), G k \A)) . (12) 

In the expression (p!2|) for the Lagrangian L, the quantity V is some appropriate domain in R 3 , 
which yet has to be specified. 

As is well known, requiring the action ( |T(]| ) to be stationary with respect to local variations 
of all the potential components A^, considered as independent quantities, yields the follwing 
field equations, 

V V {A)G» V {A) = 0,^ = 0,1, 2,3. (13) 
The "covariant gradient" V^(A) used above in Eq. (^) is a convenient notion, 

V M (A)=5 M + i^[^, ], (14) 

which will be frequently used in what follows. 

A perhaps more interesting system of equations would be obtained by coupling the gauge 
field to appropriate matter fields. Then the right hand side of Eq. (|13|) would be replaced by 
a covariantly conserved matter current. Such an addition is not absolutely essential for the 
questions pursued in this paper, and will therefore not be contemplated further here. 
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The non-Abelian Gauss law is obtained from the equations ( |T3| ) for \i = 0. Expressed in 
terms of the potential A the non-Abelian Gauss law reads as follows, 



V k (A)V k (A)A° - V k (A)A k = (15) 

Considered as an equation which determines the Lie algebra valued potential component Aq, 
for given space components A = (A 1 , A 2 , A 3 ), the equation ( (TBI) is a system of linear, elliptic 
partial differential equations with the time variable x° G [x°, x°] acting as a parameter. 

From now on it will be convenient to denote the time derivative of any quantity with a dot 
on top of that quantity, thus for example 

A k (x) = d A k (x) (16) 

The apparent equations of motion are obtained from the field equations (|13|) for \x = k = 
1,2,3. These equations read as follows, separated into two groups of first-order equations for 
convenience, 

A k = V k (A)A° - G ok (A) , G ok (A) = V e (A)G k£ (A) - tg[A , G ok (A)}. (17) 

I then return to Eq. (|15|), i.e. the Gauss' law constraint. In order to analyse the solvability 
of Eq. (pl|), one has to specify the domain V of the independent variables x = (x 1 , x 2 .x 3 ) £ R 3 
in this equation. 

In what follows, I assume that the domain V is a finite simply connected domain in R 3 with 
a smooth boundary dV. Eventually one has to contemplate the situation where V is allowed 
to grow indefinitely, but for the time being I consider only finite V. An example of a possible 
domain V would be a ball B R of radius R, centered at the origin of space coordinates, i.e. 

B R = {x| I x |< R}, (18) 

keeping in mind the possibility of the limit R — > oo at some appropriate stage in the development 
of the formalism. 

I will postpone till later the discussion of the boundary conditions which guarantee the 
uniqueness of a solution A to Eq. (0), assuming that a solution actually exists. The boundary 
conditions are related to the canonical structure of the theory, as will be seen presently. 

The solution A of Eq. (|T5|), assuming the existence of a unique solution, is, in general, an 
x-dependent functional of both the space components A and their time derivatives A, i.e. 

Ao = A {A,A} (19) 

The question is then whether the Yang-Mills system, which is originally defined by the 
action (|10D, permits a canonical structure when the potential component A is a solution to 
Gauss' law (0), i.e. when A is given in terms of A and A by the expression ([191) . This 
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question can be analysed without specifying the actual functional form of the relation (|T9"| ) in 
minute detail. 



The Lagrangian of the Yang-Mills system, when the potential component A 0} is given by 
the relation (|TP|), is obtained simply by inserting the solution (|I^) for Aq into the Lagrangian 
( 12) above. The resulting Lagrangian which will be called Lo, is then explicitly given as follows, 

L = -i^d 3 x(v fc (A)A {A,A}-A fe ,V fc (A)A {A,A}-A fc ) (20) 

-\j y d 3 x (G kl (A),G kl (A)] 



At this point it is appropriate to check whether the action principle involving the Lagrangian 
L in (|20D above reproduces the equations of motion ([17]). It is perfectly straightforward to 
verify the following result, 

8 £ Q f dx°L = - J X J dx° jVx (5A k ,\7 (A)(V k (A) A {A, A} - A fe ) -V t G M (A)) (21) 
- £j dx° d 2 a k (8 A {A, A} , V k (A)A° { A, A} - A k ) 

Here the boundary conditions for the solution Aq to Gauss' law, i.e. the system of linear 
elliptic partial differential equations (|H>D, come into play. The surface term in Eq. ([H]), if 
non-vanishing, destroys the variational principle which is supposed to lead to the equations of 
motion ( |TTD with A given by fllTf ) . It is therefore necessary to expurgate this surface term from 
the expression (|21|). This can be done by declaring that the system of linear elliptic partial 
differential equations flTBl ) for A Q has to be solved with fixed boundary values at dV, which 
results in the condition, 

<L4 {A,A}| xe9y =0 (22) 

for all admissible variations of the independent generalised coordinates A and velocities A, 
respectively. From now on it will be understood that fixed boundary values, i.e. boundary 
conditions of the Dirichlet type, will have to be used on the boundary dV in solving the Gauss' 
law (15). The appropriate Dirichlet boundary conditions will be discussed in detail in a separate 
paper 
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Under these circumstances the variational principle 



„0 

f dx°L = 0, (23) 



leads to the following equations of motion, 

V (A) {\/ k (A)A° {A, A} - A fe } - V e G ke {A) = 0, (24) 



as is evident from the relation (|2"T|). Needless to say, the equations ( |2~3| ) are nothing but the 
original equations of motion (|i~7D , albeit in a slightly disguised form, with Aq given by the 
solution (|i"9l) to Gauss' law (|T5|). 
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4 Canonical formulation of the equations of motion 



4.1 First attempt at constructing canonical momenta 

Using the Lagrangian (EDI) as a starting point, one may now attempt to construct a canonical 



description of the system defined by the equations of motion fl24|) 



The formal definition of the canonical momentum conjugate to the coordinate A a is, 

P fc V,x) = § J^ q = (v fc (A)A {A,A}) V,x) -i£(x°,x), (25) 

where the condition (B^) has been used in the calculation of the functional derivative of Lq 
in ( p5|) above. In view of the fact that Aq in Eq. ([25|) satisfies Gauss' law (15), one finds 



immediately from (|25[) that 

V fe (A)P fe (x°,x) = 0. (26) 

Now one is supposed to be able to solve Eq. (^) for the generalised velocity A% in terms of A 
and P, respectively. But this is impossible, since, Eq. (|^) can not be solved for the quantity 
T defined below, 

r(x ,x)=V fc (A)iV,x), (27) 

i.e. if one tries to derive an equation for the quantity T defined above, from Eq. (|25|), one gets 
a completely vacuous identity for this quantity, as a result of Eq. 



It would seem then, that the equations of motion (|24|) can not be written as a system of 
canonical first order equations at all. However this is actually not the case. 



It turns out, that the quantity T defined above in (|27|), can be set to zero as a consequence 
of the freedom of choosing a gauge in Yang-Mills theory. To be precise, the condition 

r(x°,x) = (28) 

has recently been shown to be a proper gauge condition, called the generalised Coulomb gauge 
condition 



In what follows, it will be assumed that the generalised Coulomb gauge condition fl2"8"| ) is in 
force. Since the generalised velocities A k then no longer are independent quantities, one cannot 
use the formula ( f25|) as it stands for the construction of the canonical momentum variables. 
An alternative procedure that leads to a canonical formalism will be given below. 



4.2 Canonical coordinates, momenta and Hamiltonian 

It is now appropriate to summarise the various conditions on the generalised coordinates A 
and genralised velocities A, respectively, that have come out of the previous analysis. 
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In the first place, it is assumed, that the generalised Coulomb gauge condition ( ffijj ) is in 
force, i.e. that 

V fc (A)iV,x) = (29) 
Then Gauss' law flT5|) simplifies, and takes the following form, 

V k {A)V k {A)A° = , x G V C R 3 (30) 

Hence i/ie solution A to Gauss' law ( |3"U[ ) is a functional of the generalised coordinates A only, 

A = A {A}. (31) 



Fixed boundary conditions, independent of A, at the boundary dV of the domain V in Eq. (|30D 
are used to define the solution ( |3~T] ) of the present form of Gauss' law fl3"U|). As a consequence 
of this requirement the following condition will always be valid, 

(L4 {A}| xe9v =0 (32) 

for any admissible variations of the generalised coordinates A. 

Finally, the Lagrangian for the system under the conditons detailed above is recorded here 
for the sake of completeness. This Lagrangian, denoted by Loo, is nothing but the Lagrangian 
of Eq. (PD|), but with the solution ([JT]) for A in stead of (|TP|). Explicitly, 



~ J v d 3 x (v k (A)A {A} - A*, V k (A)A° {A} - A k ) (33) 
d 3 x (G kl {A),G kl (A)) 



-00 



AJv 

I will now implement the gauge conditions (p9| ) as a constraint, by means of a Lie algebra, 
or rather matrix valued Lagrange multiplier field C(x), which is used to modify the Lagrangian 
L o m O) m an appropriate way, 



L' = Loo + / d 3 x(C(x), V k (A)A k ) (34) 



The Lagrange multiplier technique as such is well known in general as well as in the particular 
case of gauge systems fl2|l . 



One can now make a direct transition to a Hamiltonian using the modified Lagrangian 
above, in a manner described in the general case by Berezin ||13|. The starting point is the 



familiar definition of canonical momentum variables 7r£, 



5L' 



(V k (A)A {A}f -A a k - (V k (A)C) a , (35) 



<L4*(x°,x) 

where x° is a fixed parameter with variable x G V. The argument (x°,x) is omitted in the 



last term in (|35|) above. Similar simplifications of notation will be used freely below, whenever 



expedient. The equations (|35|) above, together with the constraint equations (|29D are now 
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supposed to be solved for the quantities A\ and C a in terms of the canonical coordinates A\ 
and momenta 7r£, respectively. 



Using Eqns. (|29|) and (0), one finds immediately from Eq. fl35|) that 

- V k (A)\7\A)C = V\A)7r k , (36) 

which, together with appropriate boundary conditions at dV, defines the quantity C as an 
^-dependent functional of A and tt, respectively, 

C = C{A,7r}(x°,x). (37) 

The boundary conditions alluded to above are very important for the consistency of the Hamil- 
tonian formalism, as will be seen presently. 

One now straightforwardly expresses the generalised velocity in terms of coordinate- and 
momentum variables, 

At = (V k (A {A}-C)T-nt. (38) 

The construction of the Hamiltonian H then proceeds in the usual way. The relation defining 
the Hamiltonian H is the following, 



H= /^ 3 x(^,i fc )-L 00 , (39) 



where the quantity A ocurring in the expressions in ([3£]) should be given in terms of canonical 
variables by the expression (^). It should be observed, that it is indeed the Lagrangian Loo 
which enters in the definition of the Hamiltonian H above, since at this stage the constraint 
(|29|) is an identity. 



By straightforward calculation one finally obtains the Hamiltonian expressed in terms of 
canonical variables from the definition (|5D|), 

H = -^d 3 x(^,7r fc ) + ^J v d 3 x (G kl (A),G k \A)) (40) 

+ j v rf 3 x (tt,, V k (A)A° {A}) + I J v d 3 x (V k (A)C, V k (A)C) , 

where the quantity C is the appropriate solution to the system of linear elliptic partial differ- 
ential equations ([$§), as discussed previously. 

By straightforward functional differentiation of the Hamiltonian H in ( ^0|) one finds, 

r tt 

^ (2;0 ' X) = ftr a *(s°,x) = ( Vk ^ Ao ^ ~ C ^ ~ ^ 
which agrees precisely with the expression (p8|) as it should. 

It is much more intricate to obtain the second half of the Hamiltonian equations, which give 
the time evolution of the canonical momentum variables % k . It turns out, that the quantity 
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C occurring in the expression (PH), which is a solution to the system of linear elliptic partial 
differential equations (|36|), must satisfy a boundary condition of the Neumann type at the 
boundary dV of the domain V in order that the functional derivative of the Hamiltonian with 
respect to the canonical coordinate variables Af, be defined. The boundary condition in question 
is the following, 

d n C\ xedv =d n A {A}\ xedv , (42) 

where d n stands for the normal derivative at the surface dV, i.e. the derivative in the direction 
of the outer normal of dV. It should be observed, that the quantity A in (f42T) is itself a 
solution to a system of partial differential equations, namely Eq. (|3T)|), for which Dirichlet type 
boundary conditions have to be used, in order to ensure the consistency of the Lagrangian 
formalism, as discussed previously. The role of the boundary conditions in the present context, 
both Dirichlet- and Neumann-, will be considerd in more detail elsewhere . 



Under the condition ( j42|) above, one finds after a lengthy but fairly straightforward calcu- 
lation, that 



y tt 

7T g V,x)^- = - ig [A {A}-C,V k (A)C + n% (43) 

V t G k \A)) -tg[C,V k (A)A {A}] a 



Needless to say, the pair of Hamiltonian equations of motion, Eqns. ((IT]) and (f43T) are 
in perfect agreement with the Lagrangian equations of motion obtained from the variation 
principle, 

5 J dx°L' = 0, (44) 

where the Lagrangian L' is given in Eq. (|34T). 

We have thus finally arrived at a canonical Hamiltonian formalism for the Yang- Mills system, 
in the case when Gauss' law is supposed to be identically satisfied. As was seen above, this 
required, among other things, that a special gauge condition, called the generalised Coulomb 



gauge condition | 10fl , was in force 



The interest in the canonical structure of Yang-Mills theory comes from the interest in 
quantising this theory, either by canonical operator methods or by functional integral methods. 
I hope to return to such questions in the near future. 

Before concluding this subsection, it is appropriate to notice, that the canonical Hamiltonian 
equations (|4l| ) and (|43T) admit a first integral, i.e. 

d (V k (A)ir k ) = 0. (45) 

The quantity Vk{A)7i k is thus a constant of motion. This is related to the fact, that the 
description of the Yang-Mills system in terms of the canonical variables A k and n% introduced 
here, is not a description in terms of a minimal number of dynamical degrees of freedom. Also 
this circumstance will be discussed in more detail elesewhere [|T1| . 
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4.3 Comparison with the Weyl gauge formalism 



For completeness I present below the Hamiltonian formulation of Yang-Mills theory in the so- 
called Weyl gauge || A = 0. This has been discussed with admirable clarity in a paper by 
Jackiw jTj, to which I refer for details. This paper also discusses the Schrodinger quantisation 
of the Weyl gauge version of Yang-Mills theory. 

The Yang-Mills Lagrangian in the case Aq = 0, which will be called Lw here, is the following, 

!/■„,: :u, 1 



- J d 3 x(A k , A k )--J d 3 x(G ke , G M ) (46) 



The Lagrangian ( |46| ) describes a theory which is not the same as Yang-Mills theory, since Gauss' 
law is absent from the field equations following from the action principle with the expression 
( [4*6]) as Lagrangian. Gauss' law takes the following form in the present case, 

V k (A)A k = (47) 

The traditional way to analyse Yang-Mills theory in the Weyl gauge, is to proceed from the 
Lagrangian fl46|), disregarding Gauss' law to begin with. 



Using the variables A^ and A^ as generalised coordinates and velocities, respectively, it is 
perfectly simple to derive a canonical Hamiltonian formalism for the system defined by the 
Lagrangian fl46|). The corresponding Hamiltonian H\y is, 



H w = -\J d 3 x(vr fc , 7c k ) + ~ J d 3 x(G M , G M ) (4* 



i.e. a simple sum of a kinetic term, depending on n only, and a potential, or interaction term, 
depending on A only. 



The Hamiltonian Hw in Eq. fl48|) should now be compared with the Hamiltonian (f40|) 
derived in the previous subsection. Disregarding for the time being the slightly different mean- 
ing of the symbols in the two Hamiltonians in question, one observes that the Hamiltonian 
( PD| ) contains interaction terms depending on the quantities A Q and C respectively, which are 
completely absent from the Weyl gauge Hamiltonian (f48|). 

A transition to quantum versions of the Hamiltonians above, by means of the fixed time 
Schrodinger quantisation rule, 

< - K = (49) 
facilitates a comparison of the two Hamiltonians. 

Interpreting the Hamiltonian ([40"D as a quantum operator H by means of the substitu- 
tion (|3S|) one is then invited to consider the following eigenvalue equation, in self-explanatory 
notation, 

HV(A) = EV(A) (50) 
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which in principle should describe a glueball spectrum. There are pitfalls associated with the 
quantisation procedure briefly described above, which have to be understood and avoided before 
considering conclusions based on equations such as fl50|) , as being ironclad. 

In the Weyl gauge case one would have, similarily, 

Hw^w(A) = EV W (A) (51) 

At this stage one simply adds the hitherto omitted Gauss' law as a condition on the 
states ^ W (A), 

1 ° 6 9f ab c A b k (x°,x) xAcr \ - W(A) = (52) 



dx k 5A a k (x°,x) " Jab fcV ' y <L4£(x°,x) 



It is now appropriate to recall the first integral fl45|) admitted by the (semiclassical) system 
described by the Hamiltonian H in (f40f ). In the Schrodinger quantisation of that system, the 



relation ( fl5D gives rise to a condition on the states ^(A) of the same type as (52), although 
not necessarily identical to that relation. Even if one might reproduce exactly a condition 
of the form (|52J) in the Schrodinger quantised theory defined by the quantum Hamiltonian 
H, by an appropriate choice of a constant of integration, one must still reckon with the fact 
that the Hamiltonian H definitely differs from the Hamiltonian Hy/ by the interaction terms 
involving the quantities A and C in the (semiclassical) expression (|40|). These interaction 
terms, appeared as a result of the demand that Gauss' law should be satisfied identically. 
The variables Aq and C in these intreraction terms had to satisfy certain very precise boundary 
conditions in order to ensure the consistency of the (semiclassical) Lagrangian and Hamiltonian 
formulation, as explained previously. Hence these interaction terms cannot be declared to 
be uneccesary or unmiportant. So, one is forced to conclude that the Weyl gauge quantum 
canonical formalism in Yang-Mills theory differs in substance from the corresponding formalism 
developed in this paper. 



5 Summary and Conclusions 

In this paper I have considered the problem of defining canonical coordinates and momenta in 
pure Yang-Mills theory, under the condition that Gauss' law is identically satisifed. Satisfying 
Gauss' law in this context means solving a system of elliptic linear partial differential equations 
for the matrix valued gauge potential component Aq, with a fixed value of the time coordinate. 
For a discussion of the canonical structure of Yang-Mills theory, when Gauss' law is identically 
satisifed, one does not need the explicit form of the solutions to the differential equations in 
question, but merely certain properties of the solutions, which are embodied in the equations 
themselves. Existence of solutions is not proved here, but assumed. The related question of 
uniqueness of solutions, which depends on boundary conditions, is only briefly touched upon. 
The requisite boundary conditions for the system of equations are not postulated a priori, but 
arise as consistency conditions related to the equations of motion. The canonical structure of 
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the theory requires the choice of a special gauge, called the generalised Coulomb gauge in this 
paper. This gauge, the existence of which is equivalent to the existence of appropriate solutions 
to a system of partial differential equations of the same type as encountered in connection with 
Gauss' law, does in principle not suffer from Gribov ambiguitites. 

The canonical variables, coordinates and momenta, as well as the Hamiltonian, are explicitly 
constructed. This construction involves again a system of partial differential equations of 
the same type as encountered previously, but now with different boundary conditions, which 
are dictated by necessary consistency conditions for the existence of canonical Hamiltonian 
equations of motion. 

Schrodinger quantisation of the theory is briefly discussed, and compared with the corre- 
sponding quantum theory in the Weyl gauge, A = 0. It is concluded that the Schrodinger 
quantised Hamiltonian formulation of Yang-Mills theory developed in this paper, differs in 
substance from the corresponding Weyl gauge formulation. 
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